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Abstract 

The effect of background dynamics of the universe on formation of large scale 
structures in the framework of Modified Newtonian Dynamics (MOND) is 
investigated. A spherical collapse model is used for modeling the formation 
of the structures. This study is done in two extreme cases: (i) assuming 
a universe with a low-density baryonic matter without any cold dark mat- 
ter and dark energy; (m) a dark energy dominated universe with baryonic 
matter, without cold dark matter. We show that for the case (ii) the struc- 
tures virialize at lower redshifts with larger radii compared to the low-density 
background universe. The dark energy slow downs the collapse of the struc- 
tures. We show that our results are compatible with recent simulations of 
the structure formation in MOND. 
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1. Introduction 



Asymptotic flat rotation curves o f spira l galaxies reveal the missin g mass 



problem in galactic scales (IBosma I . Il98ll ; Ivan Albada et al. I . Il985l ). This 



problem can be solved by cold dark matter (CDM) hypothesis, or alternative 



theories of g r avitat ion (jCaroll et al. 1 . 12004 ; lHalle et al. I . l2008t ISobouti et al. 



20091 ; IZhao 1 . 120071 ). The CDM hypothesis as a dominant paradigm together 



with the cosmologica l constant makes the standard model of cosmology 



( ISpergel et al. I . 120061 ) . Despite many efforts for direct detection of dark 



matter, it has not been observed in non-gravitational experiments. On the 
other hand, results from high resolution N-body simulations do not seem 
to be compatible with the observations on the galactic scales and also pro- 



vide incompat i ble spatial distribution o f the sub-halos (IMoore et al. I . [1999; 



Klvpin et al.l Il999l : iMetz et al. I . l2008h . 



The Modified Newtonian Dynamics (MOND) paradigm is one of the alterna- 
tive theories that was proposed by Milgrom (1983) to explain the dynamics 
of gravitational system. In MOND, the Newton's second law is modified as 
li(a/ao)a = a^ + V x H, where a^v is the Newtonian acceleration vector, a is 
the MONDian acceleration vector, a = |a| is the absolute value of MONDian 
acceleration, /i is an interpolating function which is used for transition from 
the Newtonian to the MOND regime, and an — 1.2 x 10~ 10 ms 



Dian threshold acceleration ( IBekenestein &: Milgrom 1 . 119841 ) 



" 2 is the MON- 
The value of 



the curl field H depends on the boundary conditions and the spatial mass dis- 



tribu tion, and vanishes only for some special symmetries (IBekenestein fc Milgrom 



19841 ) . This modification successfully explains the flat rotation curves of spi- 
ral galaxies at the large distances. Below a characteristic acceleration a Q , 
n(x) ~ x and in the absence of curl field H, the Newton's second law is 
approximated as a = -y/aoa^v, the so called deep MOND regime. MOND 
has rec ently been generaliz ed to a general-relativistic version which is called 



TeVeS (IBekenestein 1 . 120041 ) 



One of the challenging issues for any gravitational theory is to test its ability 
to reproduce the large scale structure formation in the universe. In MOND, 
the large scale structure formation has been studied by Sanders (1998) in a 
low-density universe, i.e. the universe without dark energy and dark matter. 
He showed that in MONDian cosmology, a patch of the universe smaller than 
the horizon size evolves with a different rate than the background. Hence 
the structures naturally can be formed through the scale dependent dynam- 
ics. In this scenario small structures form faster than larger ones (bottom-up 
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hierarchical model). Note that in this model, every point in the space can 
be assumed as a centner of collapse and the initial fluctuation of baryonic 
matter has no significant role in the formation of structures. In order to 
solve this problem, the association of MOND force can be considered as a 
peculiar acceleration in a finite- size region and the structures is formed from 
the initial baryonic fluctuation ( ISanders I . l200ll ). It has been shown that, this 
assumption leads to rapid growth of structures and the power sp ectrum of 



the structure is similar to that of obtained from ACDM universe (ISanders 



200l|). 



In the scale of clusters of galaxy, MOND reduces the discrepancy between 
dynamical and detectable mass to a factor o f 2-3 in the sca l es of clusters of 



galaxy, but does not remove it com pletely (IThe fc White I . Il988l ; ISanders 



20031 ; iPointecouteau &: Silk I . 120051 ). Sanders (2003) speculated on the re- 



maining missing mass and showed that 2-eV neutrinos which can aggregate 
on the scales of clusters of galaxy are a possible candidate to interpret the 



remaining missing mass ( ISanders 1 . 120031 ) . 



Up to present day, several attempts in the cosmological st ructure f orma- 



tion have done by means of MONDian N -body simulation (INusser 1 . 12002 ; 



Knebe et al. 1 . 12004 ; iLlinares et al. 1 . 120081 ). 



Malekjani, Rahvar and Haghi (2009, here after MRH09) showed that in 
MOND, a uniform distribution of matter finally virializes with a power law 
profile in which the virialization process takes plac e gradu ally from the center 



of structure to the outer parts (jMalekjani et al. I . l2009af ). In this work, we 



discuss the effect of background dynamics on the growth rate of the structures. 
In the case of universe dominated by cosmological constant, the spherical col- 
lapse is similar to that in CDM model, except the cosmological constant that 
changes the growth of the structure though altering the background dynam- 
ics. The main gaol of the paper is to understand the effect of background 
dynamics of the universe on the evolution of spherical over-dense structures 
in the framework of MOND. It is possible to achieve this goal even if we con- 
sider only the pure baryonic matter as a single component of matter content 
of the universe and avoid the complexity of neutrino problem in the scales 
of clusters of galaxies. Using the simple baryonic spherical collapse model 
in the context of MOND, we study the effect of generic variable dark energy 
background model on the scenario of structure formation. 
The paper is organized as follows: In section [2] we introduce the spherical col- 
lapse in MOND. In section [3] we calculate the effect of background dynamics 
on the structure formation in MOND. The paper is concluded in section HI 
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2. Spherical collapse in MOND 



The following is a brief review of MONDian cosmo logy and the spherical 



collapse in the framework of MOND (panders I . Il998l ). In the scales smaller 



than Hubble radius, the dynamics of the universe is derived from the Hubble 
equation as 

H\a) = H%{n° b a- 3 + ft°a~ 4 - (0° - l)a- 2 + (1) 

where H is the Hubble parameter at the present time, fl®, fl® and are the 
baryonic matter, radiation and total energy density parameters, respectively, 
and fide is the density parameter of dark energy. In MONDian cosmology, 
one can consider a critical radius r c that inside it the dynamics is MONDian 
(i.e., the acceleration is bellow Oo) and outside it the dynamics is Newtonian. 
The critical radius and corresponding mass scale are 

r — 2ao ]\,f _ ao 2 / \ 

' c ~ Hg| n o a -3 +2 nO a -4_ 2 Q de |; lvl c- G 'c- {*) 

The structures with mass M > M c are in the Newtonian regime, while for 
M < M c are in the MONDian regime. At the beginning of evolution, the 
dynamics of the structure is Newtonian and the density contrast grows in pro- 
portion with scale factor. Afterwards, the structure enters into the MONDian 
regime and the dynamics significantly deviates from the background expan- 
sion (see Fig. 1 of MRH09). In MONDian regime ( the scales M < M c ), 
the \x function is approximated as \i{x) = x and the dynamics of structure is 
given by 



r 

where M is the mass and r is the radius of structure. Multiplying the both 
sides of Eq.()3]) in r and takeing the integration with respect to time, we 
obtain 

v? -2(GMao)hn-, (4) 

where rj is the initial radius at the entrance time into the MONDian regime 
and Vi is the initial expansion velocity of the structure^. Changing the time 



Vi — HiTi{\ — Si) is the peculiar velocity of the structure at the entrance time to 
MOND regime, where Si is the initial density contrast at the entrance time. 
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derivative to the derivative with respect to scale factor a, Eq. (j3J) is rewritten 
as 

d i = ^f<- 2(GMo ° )lh1 ;? ( 5) 

where H(a) is given by Eq.((T]). The maximum rad ius of the structure is 
f m = r fid a ' 1 where a = v 1/^/AGMoq flSanders I . Il998h . Throughout the re- 
collapse, the global radial velocity converts to the dispersion velocity and 
eventually the structure virializes at r vir = rie 01 ^ 1 ^ 2 (for more details see 
MRH09). 



3. Background dynamics and the structure formation in MOND 



Here, we calculate the influence of the background dynamics on the MON- 
Dian structure formation for two extreme cases of low-density universe and 
dark energy dominated universe. 

3.1. Low-density background 

In the low-density background, there is no dark matter and dark energy 



i n the universe and the total density param eter is 



0.02 



(j Walker et al. I . Il99ll : ICarlberg et al.. I . Il997l ). We take three class of ob- 
jects with the masses of 1O 6 M , 1O 11 M and 1O 13 M which represent the 
globular cluster, galaxy and cluster of galaxies, respectively. By numerical 
integration of Eq. ((Sj) the evolution of structures is obtained as a function of 
redshift. The initial conditions (r^, vi), the maximum radius, r m , the virial- 



ization radius, r vir , together with corresponding redshifts 



.) are shown 



in Table (pj). The smaller structures enter the MOND regime and virialize 
earlier than the larger ones. The dynamical behavior of structure forma- 
tion in MOND is similar to the hierarchical structure formation in standard 



CDM model (IPadmanabhan I . Il993l : iMalekjani et al. I . l2009bl ). Here we see 



that the MONDian spherical collapse under a low-density background ob- 
tains the virialization at high redshift and therefore can interpret the rapid 



grow th of density perturbations. This result is also discussed in ((Sanders 



200ll ) . Sanders assumed two-field Lagrangian-based theory of MOND and ob- 
tained the rapid growth of density perturbation by solving the corresponding 
differential equation in non-linear regime. 



5 



3.2. Dark energy background 

An interesting theoretical candidate for dark energy is cosmological con 

)f state u = — lFlf Zeldovich et al. 



1967; 


Weinber . 


1989; 


Carroll . 


2001) 



However, the cosmological constant 
provides excellent fit to the SNIa and CMB data, it suffers from the coinci- 
dence and fine tuning problems at the early universe. One solution to these 
problems is assuming the dynamical dark energy models whose equation of 
state evolves with cosmic time. In these models an evolving scalar field gen- 
erates the energy and the pressure of the dark energy, and provides a positive 
accelerating universe. Here we use the variable dark en ergy model propo sed 



by Weetterich (2004) in which the equation of state is (I Wetterich 1 . 120041 ) 



where ujq is the state parameter at the present time, a is the scale factor and 
b is the binding parameter which is related to the amount of dark energy in 
the universe. However, u evolves differently for various bending parameters, 
it asymptotically converges to uq = — 1 at the present time (i.e. a = 1). In 
the variable dark energy dominated background, the density parameter of 
dark energy in Eq.([T]) can be described as 

n de (a;b,u ) = ^-ali+^o)], w^u,) = _ ^° , (7) 

where Q de = 8nGp de /3HQ and 6 = corresponds to the standard cosmo- 
logical constant model. We assume no clustering for dark energy, but it 
can influence on the growth rate of the background and consequently on the 
formation of structures. The various density parameters for variable dark 
energy model are 

n° = ng + = l-o, n° e = o.98, n b = om, n° dm = o. (8) 

Figure ([[]) shows the evolution of critical length scale function 
of redshift z, for different cosmological background models calculated from 
Eq. 02]). In the case of low-density model (black solid line), the critical 
radius increases with the scale factor. It means that the MONDian do- 
mains gets larger with time, and therefore the larger structures enter to the 



2 The equation of state is p = cup 
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MONDian regime after the smaller ones. The critical radius for the case 
of standard cosmological constant model (blue-dashed line) evolves same as 
the low-density model for z > 1, since the contribution of baryonic matter 
is dominated compare with cosmological constant in the early times. The 
cusp in blue-dashed line indicants the equality epoch of baryonic matter and 
dark energy, hence r c increases to approach infinity at this time. After the 
equality epoch, the cosmological constant is dominated and r c is smaller than 
that of low-density model. This implies that in the universe dominated by 
cosmological constant, the structure crosses r c and enters the MOND regime 
later than low-density universe. In the case of variable dark energy model 
with b = 1 (red dotted-dashed line), r c is smaller than other two models at 
any redshift. For example in the case of b = 1, the structure with the mass 
of M = 1O 11 M crosses the r c and enters the MOND regime at z enter = 68, 
while for the low-density model, this happens at z enter = 146. It should be 
noted that, because of the logarithmic term in Eq. fl6]), the critical radius for 
all models with b ^ converges to the case of b — at z *C 1. 
In order to obtain the evolution of the structure in MOND under the vari- 
able dark energy background, we integrate the Eq.(jS]) by using Eq.(jTJ. For 
different values of b, the entrance redshift, the initial radius, and the ini- 
tial expansion velocity at the entrance time for a galaxy with the mass of 
M = 10 11 M Q are shown in Table ([2]). Increasing the bending parameter 
leads the structure to enter the MOND regime at the lower redshift. Figure 
(J2J) shows the evolution of galaxy mass structure for different cosmological 
background models. The numerical results are summarized in the last four 
columns of Table (T5]). The dependence of the virialization redshift and the 
virialization radius of a galaxy mass structure to the bending parameter are 
shown in Figure (EJ). The structure virializes at lower redshift with larger 
radius for the higher values of bending parameter. The virialization in low- 
density model is equal to the virialization in variable dark energy dominated 
background with b ~ 0.4. 

3.3. Comparison with standard top-hat model 

In the framework of standard model of structure formation, the baryonic 
structures are formed within the potential well of dark matter. One of the 
simple analytical method for modeling the structure formation is the spher- 



ical collapse with top-ha t density distribution model ( IPadmanabhan 1 . 11993 



Malekjani et al. I . l2009bf ). In the standard top-hat model, the virialization 



redshift of the galaxy mass structure takes place at z vir ~ 0.5 (see Tab.(l) of 
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Ref. dMalekiani et al. I . l2009bh ). while in MOND, in the case of low-density 



background the virialaization takes place at z vir ~ 17. As we showed, the 
variable dark energy background postpones the collapse of structures. Ac- 
cording to Table fl2]), the case of b = 2 gives the virialization redshift that is 
compatible with standard top-hat model. 

3.4- Comparison with MONDian N-body simulation 

Recently, Linares et al. (2009) have presented cosmological N-body sim- 
ulation for the large scale structure evolution in the framework of MOND. 
They studied the effects of the curl-field, V x H, upon the gravitational struc- 
ture formation. They also showed that the large scale structure formation 
is faster in MOND due to strong gravity and found a correlation between 
the mass of the structure and the formation redshift in ACDM and MOND 
models. This correlation is given by 

log(M/M ) = -0.8161 \ogz vir + 11.5086, ACDM 

log(M/M ) = -0.8981 log2„ ir + 11.4685, OCBMonal (9) 

log(M/M ) = -0.6477 log z vir + 11.4036, OCBMonal2 

where OCBMond denotes to the MOND model without considering curl- 
filed in the solution of the MONDian Poisson's equation whereas OCBMond2 
is the solution with the influence of curl-field in the process of structure 
formation. These linear relations are obtained by the best-fit power laws to 
the scatter plot (see Fig. 11 of Linares et al. 2009). According to Eq. flH]), 
the evolution of structures in OCBMond2 is slightly faster than in ACDM, 
hence the high-mass objects appear to form at earlier times. There is also 
a significant difference between OCBMond and OCBMond2 with the former 
being closer to ACDM. This underlines the speed-up impact of the curl-field 
for the structure formation which is similar to the effect of dark energy. In 
Figure (T4]), we plot the mass of structure as a function of virialization redshift, 
based on Eq.(|9]). Here, we compare the results of ACDM and OCBMond2 
simulations with the results of variable dark energy background model. Since 
the curl field H is neglected in OCBMond simulation and this model can 
not justified in a fully three dimensional N-body simulation, therefore we 
do not consider OCBMond simulation in this comparison. In the case of 
OCBMond2 and ACDM the virialization redshift of galaxy mass structure 
is z v i r ~ 4 which is comparable with the virialization of galaxy in a variable 
dark energy background model with b = 0.97 (see FigJl]). Here we assume 
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that the mass of structure does not vary with time after the virialization. 
Therefore today's mass of the structure is equal with the mass of structure 
at the virialization time. This comparison can be considered as a technical 
tool for studying the background fluid of the universe. For example, here we 
see that the bending parameter b of the dark energy fluid is constrained as 
b = 0.97, using the spherical collapse in MOND. 

4. Conclusion 

In this paper we dealt with the problem of the structure formation in the 
framework of Modified Newtonian Dynamics (MOND) under the influence of 
variable dark energy dominated background model. We utilized the spher- 
ical collapse model and showed that the virialization of structures depends 
strongly on the background cosmology. We showed that the variable dark 
energy dominated background postpones the virialization of structures. Here 
the results of three background models including the low-density model, the 
standard cosmological constant model (b = 0), and the variable dark en- 
ergy model {b ^ 0) were compared. In the cosmological constant model, the 
structures virialize faster than other models (i.e., in the higher redshift). For 
variable dark energy model, increasing the bending parameter b causes the 
structure viralizes at lower redshift with larger radius. Therefore the variable 
dark energy model puts of the spherical collapse to the later times. The case 
of low-density model has an intermediate behavior such that the virialization 
redshift in this model corresponds with b = 0.4 in variable dark energy model. 
Finally, we compared the virialization of structures under the variable dark 
energy model with the recent results of MONDian N-body simulations. We 
showed that the various models of simulation are consistent with the variable 
dark energy model with different bending parameter. 

It should be noted that most of the galaxies are formed within the host en- 
vironment such as cluster of galaxies and influenced by external gravity. In 
our model, the structures assumed to be isolated system which is acceptable 
approximation when the external acceleration field is lower than the critical 
acceleration. Due to nonlinearity of Poisson's equation in MOND, the strong 
equivalence principle is violated and consequently the internal dynamics is 
affected by the external filed. Therefore MOND predicts different evolution 
for structures within the external gravity background. For studying the back- 
ground dynamics, we assumed the simple baryonic spherical collapse model 
by avoiding the complexity of neutrino problem which can be appeared in 
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the scales of clusters of galaxies. 
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Tabic 1: Numerical results for the evolution of the structure in the low-density background 
model for various mass scales in MOND. First column indicates the mass of the structure. 
Second column is the redshift of entrance time to the MOND regime. The third column 
is the initial radius of the structure at the entrance time. The forth column is the initial 
density contrast at the entrance time. Fifth column denotes the initial velocity at the 
entrance time. r m and z m are the maximum radius and corresponding redshift. z vir and 
r vir are the virialization redshift and the virialization radius of the structure. 



M 


Center 


r j [kpc] 


St x 10" 5 


Vi [km/s] 




r m [kpc] 




r vir [kpc] 


lO e M 


848 


0.033 


1.30 


14.5 


304 


0.079 


200 


0.047 


lO n M 


147 


13.4 


7.48 


315 


27.5 


48 


17.5 


28.4 


1O 13 M 


74.5 


123.7 


14.60 


1097 


8.4 


591 


4.8 


339 



Table 2: Numerical results for the evolution of galaxy mass scale (lO n M0) with the 
variable dark energy background in MOND. First column indicates the bending parameter 
of variable dark energy. The second column is the entrance redshift to the MONDian 
regime. Third column is the initial radius of the structure at the entrance time. Fourth 
column is the initial velocity of the structure at the entrance time. r m and z m are the 
maximum radius and corresponding redshift. z v i r and r„j r are the virialization redshift 
and the virialization radius of the structure. 



b 


Center 


r i [kpc] 


Vi[km/s] 




r m [kpc] 


Z v i r 


r vir [kpc] 


0.0 


146 


13.4 


285 


40 


37.5 


27 


22.7 


1.0 


68 


29.0 


401 


6 


219.2 


4 


132.9 


1.5 


45 


43.8 


442 


2.2 


510.8 


1.2 


309.8 


2.0 


36 


53.7 


475 


1.0 


916.1 


0.3 


555.6 
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6 r* — 1 — ' — 1 — ■ — 1 — 1 — i — 1 — 1 — 1 — 1 — ' — 1 — ■ — 1 — 1 — ' — 1 — 1 — 1 — ' — r 




-6' 1 

2 1 0-1-2-3 

log z 

Figure 1: Log-log plot of the MONDian critical radius, r c , as a function of redshift 
for various cosmological background models as described in legend. The cusp on 
the blue dashed line is due to the equality epoch between baryonic matter and 
dark energy, hence r c increases to infinity at this time. 
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Figure 2: The evolution of a structure with the mass of (10 11 M Q ) from the en- 
trance redshift, z en t er until the virialization stage for different values of bending 
parameter. Increasing the bending parameter causes the slower collapse. The case 
of b = refers to standard cosmological constant background. For comparison, 
the evolution of the structure in the low-density background is potted (black solid 
line) . 
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Figure 3: The dependency of the virialization redshift (upper panel) and the virial- 
ization radius (lower panel) of 10 11 M Q structure as a function of bending parameter 
in variable dark energy dominated background model. By increasing the bending 
parameter the structure virializes at lower redshift with larger radius. The hor- 
izontal dashed line indicates the virialization in low-density background model. 
The virialization in low-density model is equal with the virialization in variable 
dark energy background with b ~ 0.4. 
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Figure 4: The mass of structure in logarithmic scale as a function of virialization 
redshift, based on Eq. Q. The horizontal dotted-dashed line represents the galaxy 
mass structure, the vertical dashed line indicates the virialization redshift of galaxy 
mass structure in OCBMond2, ACDM. The virialization of a galaxy in ACDM and 
OCBMond2 is equal with that of obtained in a variable dark energy model with 
b = 0.97. 
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